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Abstract 

Within the standard model, we calculate the radiative B — ► K*j decay rate based on a 
Bethe-Salpeter description for the meson wave functions and the hadronic matrix elements. 
With a reasonable choice of parameters the branching ratio BR(£? — ► K*"f) is found to be 
(3.8 - 4.6) x 1(T 5 , which is in agreement with the CLEO experimental data. We also find 
with m b = 5.12 GeV the ratio R = T(B — ► K*i)[T{b — ► sj) = (10 - 12)%, which can 
be slightly larger if a smaller m& is chosen. In this approach, the light degrees of freedom 
in mesons are treated as light constituent quarks with relativistic kinematics, and the form 
factors in the decay amplitude are essentially determined by the relativistic kenematics and 
the overlap of wave functions of the initial and final state mesons. Due to the large recoil 
momentum of the K* meson at the B meson rest frame, the form factors are sensitive to the 
overlap integral of the meson wave functions, which are determined dynamically by a QCD- 
motivated inter-quark potential. Relativistic effects on the meson wave functions mainly due 
to the Breit-Fermi interactions are found to be significant in determining the decay rate. 



I. INTRODUCTION 

The interest in studying rare B decays lies with the fact that these decays, induced by the 
flavor changing b — > sj neutral currents, are controlled by the one-loop electromagnetic penguin 
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diagrams. They play important roles in testing loop effects in the standard model SU(3)c X 
SU(2)l x U(l) and in searching for the so-called "New Physics" because they involve many 
important standard model parameters such as the top quark mass m t and the Cabbibo-Kobayashi- 
Maskawa matrix elements \V t b\ and \V ts \. 

It has been shown that at the quark level the inclusive decay b — > s'j has a sizable enhancement 
due to the QCD corrections0. However, it is unfortunate that we can only observe the exclusive 
channels, such as B — > K*j, which are hadron transitions and therefore plagued by uncertainties 
in determining the weak hadronic form factors. Recently the CLEO Collaboration reported the 
result @ 

BR(B — ► K*i) = (4.5 ± 1.5 ± 0.9) x 10" 5 , (1) 

which makes it possible to test various models including the standard model and models involving 
new physics, provided that the estimate of hadronic matrix elements is under control. 

By now, there have been manv methods to calculate the hadronic matrix elements, such as the 
nonrelativistic quark model @ @ @, the HQET( Heavy Quark Effective Theory) where the strange 
quark s is assumed to be heavyil , the HQET combined with chiral symmetries where the s is con- 
sidered lig htfl , and the QCD sum rulesS. The predicted value of R = T(B — > K*^f) /T(b — > 57) 
is very different among various models, ranging from 4.5% B to 40% Q. Theoretically, the diffi- 
culty is mainly due to the large recoil momentum of the K* meson in this B decay process. 

As an attempt to tackle this problem and, in particular, to incorporate the relativistic effects 
of the underlying dynamics, in this paper we will use the Bethe-Salpeter description for hadronic 
form factors and meson wave functions^ to calculate the hadronic matrix elements involved 
in this rare B decay. In this approach, the light degrees of freedom in mesons are treated as 
light constituent quarks with relativistic kinematics, and the form factors in the decay amplitude 
are essentially determined by the relativistic kenematics and the overlap of wave functions of the 
initial and final state mesons. The form factors may then be related to the inter-quark dynamics 
through the meson wave functions. This method might allow us to calculate the form factors in a 
wide range of values of the squared momentum transfer, q 2 , provided that a good understanding 
about the covariant decay amplitude and BS wave functions is achieved. The remainder of this 
paper is organized as follows. In Sec. II we review the method we use. In Sec. Ill we apply our 
method to calculate the rate of the radiative decay B — > K*^ and compare it with those given 
by the nonrelativistic Schrodinger equation and the naive scaling law . Conclusions are given in 
the last section. 

II. FORMALISM 

The standard model SU(3)c x SU(2)l x U(l) has gained great successes and it has now been 
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thought as the fundamental theory to describe the weak-electromagnetic and strong interactions. 
However, there is one problem which frequently obscures theoretical predictions in the standard 
model. That is the so-called nonperturbative(long distance) effects that cannot be delt with 
effectively at present. 

Nonperturbative effects play significant roles in the weak decays of the heavy flavor mesons. 
Dealing with the hadronic matrix elements in these decays, one assumes that perturbative short 
distance(~ ^U) effects and nonperturbative long distance (~ Aq GD ) e ^ ec ^ s can ^ e treated sepa- 
rately. Perturbative effects can be calculated using the well known short distance techniques of 
QCD as modifications to the weak Halmitonian. Nonperturbative effects(including the exchange 
of soft gluons, the creation of quark- ant iquark pairs from vacuum and the final state interactions) 
are absorbed into the initial and final hadronic states. When dealing with the hadronic matrix 
elements one must have a good mastery of the hadron wave functions connected with the QCD 
nonperturbative dynamics. Unfortunately, at present nonperturbative effects are still difficult to 
be calculated from first principles. Therefore one often relies on various phenomenological models 
such as the Bauer-Stech- Wirbel(BSW)modelEJ, the Isgur-Scora-Grinstein-Wise nonrelativistic 
quark model@, etc. 

In quantum field theory, a basic description for the bound states is the Bethe-Salpeter equation!! 
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Define the Bethe-Salpeter wave function of the bound state | P) with an overall momentum 
P of a quark if)(xi) and an antiquark ip(x 2 ) 

X (x 1 ,x 2 ) = (0\Til){x 1 )$(x 2 )\P), (2) 

where T represents time-order product, and transform it into the momentum space 

X p(q) = e~ iPX ( d 4 xe-^ X {x ll x 2 ). (3) 



Here we use the standard center of mass and relative variables 

X = rjtxx + 772X2, x = xi — x 2 , (4) 
where T]i = j^j^a {i — 1,2). Then in momentum space the bound state BS equation reads 

{A - m l ) XP (q)^2 + m *) = ^J d ' kG (^ ? - k)xp(k), (5) 

where pi and p 2 represent the momenta of quark and antiquark respectively, G(P, q — k) is the 
interaction kernel which dominates the inter-quark dynamics. According to Eq. ([|) we have 

Pi = rjxP + q, p 2 = rj 2 P - q. (6) 
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Note that in Eq. ([5]) rti\ and rri2 represent the effective constituent quark masses so that we 
could use the effective free propagators of quarks instead of the full propagators. This is an 
important approximation and simplification for light quarks. Furthermore, because of the lack 
of a fundamental description for the nonperturbative QCD dynamics, we have to make some 
approximations for the interaction kernel of quarks. 

i) To solve Eq. fl5|) one must have a good command of the potential between two quarks. 
However, the reliable information about the potential only comes from the lattice QCD result, 
which shows that the potential for a heavy quark-antiquark pair QQ in the static limit is well 
described by a long-ranged linear confining potential ( Lorentz scalar Vs ) and a short-ranged one 
gluon exchange potential ( Lorentz vector Vy ), i.e,&3 

_i _i, 4 a (r) 

V s (r) = \r, V v (r) = ---±J-, (7) 

6 r 

The lattice QCD result for the QQ potential is strongly supported by the heavy quarkonium 
spectroscopy including both spin-independent and spin-dependent effects. In the next section we 
will employ the potential below regardless of whether the quarks are heavy or not 

V(r) = V s (r)+ llx ®YVv(r), 

Vs(r) = Xr- 

ar 

V v (r) = -l^ e — , (8) 
6 r 

where the introduction of the factor e~ ar is to avoid the infrared(IR) divergence and also to 
incorporate the color screening effects of the dynamical light quark pairs on the "quenched" QQ 
potential ®. It is clear that when ar <C 1 the potentials given in (^) become identical with that 
given in (^). In momentum space 

G(p) = Gs(P)+i»®rG v (p), 
Gs(P) = ~-5 3 (P) + , 

Ov0) - A (9) 
where a s (P) is the well known running coupling constant and is assumed to become a constant of 
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0(1) as V -> 

12?r 1 

a s (P) = ^7 — ■ ( 10 ) 



27 



ln(a 



A 



QCD 

The constants A, a, a, and Aqcd are the parameters that characterize the potential. In the 
computation in the next section we will use 

A = 0.183GeV 2 , a = OMGeV, a = e = 2.7183, A QCD = 0.15GeV, (11) 

and will discuss the sensitivity of the results to the values of parameters later. 

ii) In solving Eq. fl5|), in order to avoid the notorious problem due to the excitation of the 
relative time variable we have to employ the "instantaneous approximation" . Meanwhile, we will 
neglect the negative energy projectors in the quark propagators because in general the negative 
energy projector only contributes to quantities of higher orders in — , where tuq represents the 
mass of the heavy quark. 

We write therefore in the "instantaneous approximation" 

{i>\ ~ mi)xp(g)(*fe + m 2 ) = ^ / d^kGiP, q - k) X p(k), (12) 

where G(P, Q — k) represents the "instantaneous" part of the potential G(P, q — k). This suggests 
the derivation of an equation for the three dimensional BS wave function 

= jdq\ P {q\q) (13) 

by dividing both sides of Eq. (|T2"| ) by the propagators of two quarks, then integrating over k° and 
neglecting the negative energy projectors 



where 



^p^) = po_ El _E 2 Al ^° I d3kG ( p > q ~ fc)*-(fc)7°A 2 -, (1 I) 



A+ = 2£t( E i + 7° 7 ■ Pi +m l7 °), 

A 2 _ = t^t(E2 - 7° 7 ■ P2 -m 27 °), (15) 
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are the remaining positive energy projectors of the quark and antiquark respectively. Here, E\ = 
m\+ pi , E 2 — \J T>2 ■ From Eq. ([TJ]) it is easy to see that 

A^(?) = 

$ ? (g)A 2 _ = *f(q). (16) 

Considering the constraint of Eqs. (|16|), and the requirement of space reflection of 7°<3>^(q , )7° = 
— $ ->■(— Q ) for a negative parity meson, and the constraint of the general form of the meson wave 
funstion in the rest frame (see Eq.(19) below), then the wave function $^(?) in a moving frame 
can be written as follows'0 

d>°_T(g) = Al 7 °(l + A 75 7°A 2 M?) = ^^(1 + —h / 2 ~ m2 y>-(g), 

$L-(g) = AW(1 + Z) ^ 7 o A 2 /.(g) = ^ + mi (i + Z) / 2 ~ m2 /-(g), (17) 
p v ; +' v M y ->p^ 2E1 M 1 ^ 2E 2 p 

where $^ (<?) and $L (g) are the three dimensional BS wave functions of the 0~ meson and 1~ (S- 
wave) meson respectively. P M and M are the 4-momentum and mass of the meson. ^ = 7 M e M , e M is 
the polarization vector of 1~ meson. <p^(Q) and f-(Q) are scalar functions of P and ? in general. 
Note that in Eq.(17) the appearance of the positive energy projectors for quark and antiquark in 
the meson wave functions at the moving frame is an immediate consequence of Eq.(16). It is easy 
to see that if taking the heavy quark limit mi — > oo, then p^ — > P M , Eq. (|17"D becomes 



p v u 

*£W = ^(l+jO ^°A 2 _/«(?), (18) 

where f M = and y?„ = / v , which is due to vanishing color-magnetic force in the heavy quark 
limit. This indicates that in the heavy quark limit the BS wave functions respect the flavor-spin 
symmetry, and the light degrees of freedom are described by (p v , the wave function of the light 
constituent quark, which is to be determined dynamically by the BS equation. In the rest frame 
of the meson (P = 0) 

$ ? (g)=A^ 7 (l + 7 ) 757 °A 2 _^), 

= Ah°(l + 7°) *y°A*/(?). (19) 



It is easy to showli that Eq. (|19|) is the most genernal form for the 0~ and 1" (S-wave) q{q2 meson 
wave functions at the rest frame (e.g. for the 0~ meson wave function there are four independent 
scalar functions but with the constraint of Eq.(16) those scalar functions can be reduced to one 
and expressed exactly as Eq.(19)), and Eq. (|17|) may be obtained by boosting the spinors from 
the rest frame to the moving frame. 

Substituting Eq. (|19"D into Eq. (0), one derives the equations for <p(Q) and /(<?) in the meson 
rest frame 



Mtp x (q) = {E 1 + E 2 )i Pl {q) 

' d 3 k{G s (q -k)-AG v (q -k))fi(k) 



where 



-±2 

E X E 2 + m 1 m 2 + Q 



±E X E 2 



AE\E 2 J E\ 4~ E 2 

Ei + E 2 f 3 ^ ^ N mi+m 2 



+ T^ir / d3A;(G5( ^ - ^ + 2Gy( ^ - * ))( * ■ ^^^m/ 1 ^' (20) 

= A El E 2 ( mi+m2 ) (21) 



Mf x (q) = (E 1 + E 2 )f 1 (q) 

-— l — I d 3 k(G s (q -k)~ 2G v (q - k)){E 1 m 2 + E 2 m x )h(k) 
1 ±hj\hj 2 J 

Ei + E 2 r 3 - E^ + E 2 m x -\ 
-^E-J d kGs{Q ~ k) m 1+ m 2 fM 

+ El E 2 - mim2 + q 2 r d , k{Gs{q AGv{q _ t))Cq t)hCk) 

AE X E 2 q J 

_ E 1 m 2 -E 2 n ll r ^ ? _ ^ ^ E^ ^ 

4E 1 E 2 q J m 2 + mi 

E 1 + E 2 -m 2 - mi r ? _ - ? _ ^ . _ 1 . fS) 



2E 1 E 2 q J E 1 + E 2 + m 1 + m 2 
+ f ^4tf-i) 2 r , - 1 /:!/,•!. (22) 



£i-E 2 9 .Bi + ^ 2 + mi + m 2 ' 
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where 

/i00 = J^"^ (23) 

In the nonrelativistic limit for both quark and antiquark , Eq. (|20|) and fl2~2|) can be expanded 

in terms of 9 /mi and 9 /m 2 , and they are identical with the Schrodinger equation to the 
zeroth order, and with the Breit equation to the first order. If the antiquark is light and becomes 
relativistic then Eq. ( pop and (^) will include the higher order relativistic corrections. These 
equations will be solved numerically. 

In the BS description the transition matrix element for | P, M) — >\ P', M') is given by 

(P', M> | QTQ | P, M) = (2n)H J d 4 p 2 Tr {xP'(q')T X p(q)(^ + m 2 )} , (24) 

where Xp'W) = ~l°Xp'{q')l° ■ Eq. (|24"D can be reduced into a simpler form expressed in terms 
of the three dimensional BS wave functions, on condition that the negative energy projectors in 
quark propagators are neglected and the kernel is independent of the relative time variable 

(P',M' | QTQ \P,M) = (2tt) 3 / d 3 p 2 Tr |$L(g')7°r$ ? (9)| , (25) 

where the quark may change its flavor while the antiquark remains a spectator (see Fig. 1). Then 
the normalization of the wave function Q^(q) reads 



(2tt) 3 J d z qTr |$L(g)$-(g)| = 2E = 2\j M 2 + "p (26) 

Based on the formalism above with the solutions for the BS wave functions , we can calculate 
the form factors in various processes. We can also calculate the decay constants and the mass 
differences between 0~ and 1~ mesons, etc. Along these lines a rather extensive investigation has 
been made and the results are found in agreement with the experiments®. In the next section 
we'll apply this formalism to compute the decay rate of the radiative B — > K*^ decay. 

III. The Decay B — ► K*-f Within the Standard Model 

Within the standard model, the inclusive b — > s'j decay is governed by the electromagnetic 
penguin operator, for m s nib ^ 

H eff = Cm b e"sa^q v {l + 75 )6, (27) 



S 



where e M and q^ are the polarization vector and momentum of the photon. The constant C includes 
the QCD corrections'!!) and the dependence upon the CKM matrix elements and the heavy quark 



masses 



where 



C 



G F e 



mr 



m 



2 II 
W 



Fo(x) = r 23 I F 2 (x) + 



116 
~2T 



1/111 -\ 1/28 . 

-(r23 — 1) -| (r23 — 1 

5 V ; 14 V 1 



(28) 
(29) 



with r = 

a s (m W ) 



and F 2 (x) given by 
F 2 (x) 



x 



x-1) 



7 3x 2 - 2x 



3 12 12 2(x - 1) 



In a; 



(30) 



The function F 2 (x) depends weakly upon the top quark mass with a value increasing from 0.55 
to 0.68 in the range 90GeV <m t < 210GeV. 
The amplitude for B(P) — > K*(k,r]) , y(q,e) is 



A (B — > iT 7 ) = Cm b {K*{k, V ) \ sa^{l + 7 s)& I B(P))e?<t 



(31) 



where rj and k are the polarization vector and momentum of the K* meson and P is the momentum 
of the B meson. The hadronic matrix element involved in this process can be expressed in term 
of its Lorentz structures as follows(or equivalently with form factors fi, f 2 and /s)® 

(K*(k, rj) | sa^b | B{P)) = e^{Arf a P p + Bifjzp + Crf ■ PP a k„), (32) 

where the form factors A, B and C are functions of q 2 = (P — k) 2 . Using the Dirac matrix identity 

^75 = ~ie^a aP , (33) 
the matrix element of the current sa^^b is given by the same form factors 



(K*(k, rj) | sa^^b | B(P)) = i[ A{r, ni P v - rf^P 1 *) + B^k" - ri* u k") + Crf ■ P{P»k v - P u k») ]. 

(34) 

For the real photon these form factors are evaluated at q 2 = and the form factor C does not 
contribute to this transition. The decay width calculation is straightforward, 



r(B — ► K*i) 



\C\ 2 (m 2 B — m 2 Ki ,) 3 m 2 



47rm 3 



\A(0)+B(0)f 



(35) 



B 
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We will calculate the form factors A and B at q 2 = in our BS formalism. It is convenient to 
do the computation in the rest frame of B meson (P= 0) in which the wave functions of B~(bu) 
and K*~(su) can be expressed as follows 

$b(?) = Tyrri (1+7 )75^7^ ¥>(?), 

♦.•(« ) (36) 

where 9 and ? are the internal relative momentum of B meson and K* meson respectively: 
Q=-Pu,Q= k-Pu- Substituting Eqs. © into Eq. (H) we get 



(K*(k, rj) | sa^b | B{P)) = {2ixf J d 3 p u Tr{yj(l + 



9„ + m a fl 



M K * ' 2E S 



7 <V7s 



9F (1 + 7 MPu)/( fc -Pu), (37) 

where the u is a spectator in the transition, and 



E s = \lm 2 s +p s 2 = \lm 2 s + (k - Puf 



E b = \j m 2 b + p b 2 = \j ml+ Pu , E u = \Jml+p u 2 . (38) 
Here we have employed the identity 

Pu ~ m u ^ i> u -m u = j> u -m u 
2E U 2E U 2E U 

Considering the form factors are independent of the polarization of K* and Lorentz indices, we 
may choose some specific polarization(e.g. a transverse one) and Lorentz indices to do calculation. 
From Eq. (^) We obtain 

D ( 27F ) 3 f J3 rr (jifi , H \^s + m s o i i A + m b 



Tftf-rftfJ L " v Mr-* 2E s 111 10 2E, 

(1 + 7%^^}*(a)/(— ^ * - A), (40) 
2E U m s + m u 

A=M (2,f[£p u Tr {^(l + ^) A + m - 7^7 5 A + mfe (l + 7°)75 



"Mp»)/( r~ ifc - ft,) - B J Ml* + P ]. (41) 



2E U m s + m u 
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Apart from two approximations i.e. neglecting the dependence of the kernel on the relative 
time and neglecting the contribution of higher order negative energy projectors, the expressions 
( fiOD and ( f4"T| ) are rather genernal and model independent. To calculate the values of A and B we 



need to know the scalar wave functions <p and /, which depend on the dynamical model to be used. 
Here we will use the potential model described in (^) — ([TT| ) to solve for ip> and /. Substituting ip 
and / obtained by solving Eqs. ©,(0) and Eqs. (||),(||) into Eq. (|0]) and Eq. (0) we get 

|A(0) + 5(0)| = 0.62 (42) 

and 

55(5 — ► K*i) = 3.8 x HT 5 , (43) 

where we have used \V ts \ = 0.042, m t = 150GeV, \V tb \ ~ 1, m b = h.\2GeV, m s = 0.55GeV, 
m u = 0.33GeV and t b ^ 1.3ps. Our result is close to those of Ref.0 (|A(0) + 5(0) | = 0.53) and 
Ref.|fl6f (|yl(0) + 5(0) | = 0.46 in our notation). Taking into account the experimental uncertain- 
ties of | Vfs | (ranging from 0.030 to 0.0540), we find the branching ratio 55(5 — > K*^f) ranges 
from 2 x 10~ 5 to 6 x 10~ 5 which agrees with the CLEO data given in (|I|) within errors. The ratio 
of the exclusive rate to inclusive rate, R, is expressed in term of A and B as follows 



Substituting Eq. (g2|) into Eq. (||) we get 

R = 10%, (45) 

In comparision we have solved the nonrelativistic Schrodinger equations for the scalar wave 
functions (p and / in the nonrelativistic limit using the same potential as that in the BS equations 
i.e. Eqs. (EDj) — ( E51 ) which are automatically reduced to the nonrelativistic Schrodinger equa- 
tions to the lowest order in q /m?(i = 1,2). We find the value of \A(0) + 5(0)| is about 0.32, 
significantly smaller than that given by Eq. fl4"2"|) . 

Meanwhile, we have also assumed the scalar wave functions are the Gaussian as in some 
nonrelativistic quark model Hfltl 

(p(S) oc e-^ /a \ f(q) oc l h \ (46) 

There the parameters a and 6 were obtained by the variational method in the nonrelativistic 
Schrodinger equation @ . Here we will use a simpler method to estimate their values. The param- 
eters a and b are connected with the mean value of the internal momentum squared of the 5 and 
K* mesons 

n 4 ^2 4 ^2 

a 2 = -(q) B , b 2 = -{q) K *. (47) 
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Using the "virial theorem" 

and extrapolating the naive scaling law obtained from the study of the heavy quarkonium: V(r) > 



Clnr where C = 0.73GeV to the B and K* mesons, we get 

( t) = txC, (49) 
where /i represents the reduced mass of the meson system. Therefore, for B meson and K* meson 

a 2 = 0.30GeU 2 , b 2 = 0.20GW 2 , (50) 

which are consistent with those obtained in the nonrelativistic quark models! i i. Using 
Eqs. (H), fl5Q|) and Eqs. ©,([41]), we find the value of \A(0) + B(0)\ is about 0.31, almost in 



coincidence with our result based on the nonrelativistivic Schrodinger equation solutions. 

From the results we obtained above we see that the decay rates given by the relativistic BS wave 
functions distinguish them significantly from those given by the nonrelativistic Schrodinger wave 
functions. This reflects the dynamical differences between these two descriptions and suggests 
that the relativistic effects must be considered when we deal with the systems containing light 
quarks. The scalar functions of different desciptions are compared in Fig. 2 and Fig. 3 . From Fig. 
2 and Fig. 3 we see that the wave function of B meson in the BS decription is "fatter" than that 
in the nonrelativistic Schrodinger decription and that the 0~ meson wave function has a longer 
"tail" than the 1~ meson. These are relativistic effects which are mainly due to the well known 
Breit-Fermi interactions, including both spin- dependent and spin- independent terms. In fact, in 
the 0~ channel there is a very short-ranged (5 function like) spin-spin force between quarks, which 
is attractive and lowers the energy level of the 0~ state and pulls the quarks towards the origin. 
Consequently, in momentum space the 0~ meson wave function becomes "fatter" and has a longer 
"tail". Whereas for the 1~ mesons the spin-spin force is repulsive and weak (three times weaker 
than that for the corresponding 0~ mesons), and is compensated and even overwhelmed by other 
attractive spin-independent relativistic corrections. These dynamical ingredients are contained 
in Eq. (BO) for the 0~ meson with Eq. ( |22] ) for the 1~ meson, and can be explicitly seen by the 

^ 2 9 ^ 2 

nonrelativistic reduction of these equations in terms of Q /m\ and 1 /m-z- 

In line with the observation made in Ref||, we find that our results depend rather strongly 
upon the wave functions. When we change the "characteristic momentum" parameters a and b in 
Eq. ((4(]), we find that the value of \A(0) + -B(0)| changes rapidly. This is because the momentum 
transferred to K* is very large and K* is so far away from the zero recoil limit that the decay rate 
essentially depends upon the overlap of the wave functions of the initial and final meson states, 
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and in this kinematic region the overlap becomes particularly sensitive to the broadness of the 
wave functions in the momentum space. The more broad the wave functions are, the larger the 
overlap and hence the decay rate. 

The numerical values of branching ratio (43) and ratio R (45) are obtained by using the 
potential parameters (11) and quark mass parameters given after (43). It is significant to examine 
the sensitivity of the results to those parameters. In doing so, we first use the same quark 
parameters as before but change the potential parameters. We find that the value of |A(0) + -B(0)| 
is insensitive to a and a (see (11)) (note that the screening effect of a is mainly on higher excited 
states rather than the ground state mesons, and that the strength of the running Coulomb force at 
large distances, which is associated with a, also has little effects on the ground state heavy mesons). 
On the other hand, \A(0) + 5(0) | is increased as the string tension A and Aqcd increase. This is 
because, larger A and Aqcd result in stronger attractive inter-quark forces, and therefore make 
the meson wave funstions more compact in coordinate space and more broad in momentum space, 
and hence increase the overlap integral of meson wave functions and the value of \A(0) + B(0) \. 
With a popular choice of A = (0.18 - 0.20) GeV 2 and A QCD = (0.15 - 0.20) GeV and other 
parameters unchanged, we find instead of (43) and (45) 

BR(B — ► K*j) = (3.8 - 4.6) x 10~ 5 , (51) 

and 

T(B — > K*j) , 

R ^ r ( ^., 7 ) U (io-i2)%. (52) 

Next, we use the same potential parameters as (11), but change the quark mass parameters. We 
find that \A(0) + B(0)\ is insensitive to the b quark mass. This is because, as naively shown in 
(49), the internal momenta of quarks in a meson are mainly determined by the reduced mass. 
With m u = 0.33 GeV, m b = (4.70 - 5.12) GeV, the reduced mass // = (0.308 - 0.310) GeV is 
almost unchanged, and hence \A(Q) + B(0)\ remains rather stable. However, from (44) it is clear 
that the ratio R is sensitive to mj. With mj going down to 4.70 GeV from 5.12 GeV, R will be 
increased by about 30% entirely due to a smaller in the mass ratios in R, while R is decreased 
by less than 10% due to a slightly smaller value of \A(0) + B(0)\. It is obvious that a larger value 
of R will favour a smaller value of m^. 

Very recently the CLEO Collaboration has reported the first result of inclusive b — > 37 decay 



branching ratio 



BR{b — ► ai) = (2.32 ± 0.51 ± 0.29 ± 0.22) x 10" 4 , (53) 

T(B — ► K^) , . 

i? = -^— = 0.19 ±0.08. 54 

r(6 — > s 7 v ; 
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Our predicted value (52) for R is smaller than the experimental value (54). However, as discussed 
above, if we use a smaller value for m&, the calculated value for R can be increased by, say, about 
(10 — 20)%. Nevertheless, it seems to be difficult for our model to reach the large central value of 
R = 0.19 measured by CLEO. Further investigations are still needed. 



IV. SUMMARY and CONCLUSION 

In this paper we solve the BS equations for mesons and employ the covariant form for the 
wave functions and the transition matrix elements to calculate the form factors involved in the 
radiative B — ► K*^ decay. In principle, we can calculate the form factors at any values of 
the squared momentum transfer. In practice, we have made two approximations, i.e., neglecting 
the dependence of the kernel on the relative time and neglecting the contribution of the higher 
order negative energy projectors. We obtain a rather genernal form for the form factors ( flUD 
and (0). These expressions are rather model independent ( apart from the two approximations 
mentioned above). We then use a QCD- motivated one gluon exchange plus linear confinement 
potential as the kernel to solve the BS equation. With a rather popular choice for the potential 
parameters and quark masses, we get the scalar wave functions ip and /, and then calculate the 
form factors A and B. For \V ts \ = 0.042, \V t b\ = l,m t = 150GeV we find that the branching ratio 
BR(B — ► K*7) = (3.8 - 4.6) x 10~ 5 , and R = (10 - 12)%. R may take a slightly larger value of 
(12 — 14)% if, instead of choosing m h = 5.12 GeV, a smaller b quark mass say m b = (4.7 — 4.9) GeV 
is used . 

Because B — ► K*j is a large recoil process, the relativistic effects are important. There are 



two sources of relativistic effects. One is from relativistic kinematics, which may be seen in fl40D and 
(p|) where the form factors are expressed in terms of Dirac spinors, and where E u 3> m u , E s ^> m s . 

In (0) and fl4T|), due to the large recoil momentum of the K* meson | k \ — ( - Ma 2 ^ g * ^ — ^f 2 -, 
the overlap integral of the wave functions of B and K* becomes much smaller than that in the 
zero recoil limit. In connection with this , there is another source of relativistic effects, i.e., the 
dynamical effect on the meson wave functions. With the large recoil momentum the overlap is 
small and therefore is particularly sensitive to the meson wave functions, which are determined 
by inter-quark forces. To see this we have used three forms for f and / in ( |40|) and (f4T|) : (a) the 
solutions of BS equation; (b) the solutions of zeroth order Schrodinger equation with the same 
inter-quark potential; (c) the Gaussian wave function. We find that the three results are quite 
different and the relativistic effects on the wave functions (p and / are indeed important. With 
(b) and (c) the obtained decay rates are similar but smaller by more than a factor of 2 than 
with (a). This is because, with BS equation the Breit-Fermi interactions induced by relativistic 
motion will broaden the wave functions in momentum space. The important effect is that the B 
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meson wave function is broadened by the color magnetic i.e. the hyperfine spin-spin force and 
other spin-independent terms induced by one gluon exchange, which also gives the B* — B mass 
splitting with a right size@. With the broadened wave functions the overlap integral is increased 
and a larger decay rate is achieved. In this connection, the hadronic matrix element involved in 
the decay B — > K*^ may indeed provide a test of some important ingredients in the inter-quark 
dynamics. The effects of Breit-Fermi Hamiltonian on the wave functions have also been shown in 
charmonium decays®. The observed suppressions for the electric dipole transitions ip' — > •yxj 
and xj — ► 7^/V ; (J = 0,l,2), as well as l P\ — > jrjc are probably due to these effects, because 
the broadened wave functions in momentum space (hence narrowed in coordinate space) by the 
Breit-Fermi interactions will reduce the dipole transition overlap integrals, and therefore the rates. 

Of course, there are theorectical uncertainties in our approach, such as the neglect of the 
retardation effects and the contribution of the negative energy projectors, the lack of knowledge for 
the correction to the static inter-quark potential due to the light quark motion, and the gluon(hard 
and soft) exchange nonspectator effects that are difficult to compute in the quark models at 
present. Hopefully, our results can serve as an useful estimate of this decay. Definite conclusions 
strongly depend upon the reduction of these uncertainties in the theoretical computations which, 
we hope, will be improved in the future. 
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Note added. After this work was completed and submitted for publication the CLEO Collab- 
oration reported a result of the inclusive b — *> s'j decay rate (Ref. [19]). We have therefore included 
a brief dissussion for it in our calculations. 
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Figure Captions 

Figure 1: Diagram for the meson transition | P, M) — >\ P',M'). Here the quark changes its 
flavor and momentum via the vertex T, while the antiquark remains a spectator. 

Figure 2: The B meson wave function <p(P), defined in Eq. (plj) and normalized by Eq. ([26|). 
The solid line represents the BS equation solution by solving Eqs. (^t]) and (^T|) while the dashed 
line represents the Schrodinger equation solution by solving Eqs. ( ^0|) and ([H]) but only keeping 
the lowest order terms in the nonrelativistic expression. 

Figure 3: The K* meson wave function f(P), defined in Eq. (|19|) and normalized by Eq. (Effi). 
The solid line represents the BS equation solution by solving Eqs. (|^) and ([2^) while the dashed 
line represents the Schrodinger equation solution by solving Eqs. ( p2"l) and (^) but only keeping 
the lowest order terms in the nonrelativistic expression. 
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